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SUMMARY 



The flow of a compressil)ie a sphere fixed 

in a uniform stream is calculated to the third'order of 
appro-ximation "by means of the .Janzen-Raylei^?h method. 
The velocity and the pressui'e dlstri'but.ions over the sur« 
faqe of the sphere are computed and the terms involving 
the fourth power of the Mach numher, neglected in Rayleiajh 
calculation, are shown to he of considerahle importance as 
.the local velocity of sound is approached on the sphere. 
The critical Mach numher, that is» the value of the Mach ■ 
number at which the maximum velocity of the fluid past, the 
sphere is just equal to the local velocity of sound, is 
calculated for "both the second and the third approxima-^ 
tions and is found to he respectively , M^j. = 0,587 and 
0.575, 



The ^.rrotational flow of a compressi fluid past a 
circular cylinder and a sphere was first calculated hy 
Janzen (reference l) and hy Raylei^h (reference 2). Their 
method consisted in ohtainin^ va correction term to the in- 
corapressihle-fluid solution, "but the results were limited 
to the terms involving only the ^^quare of fhe Mach number* 
Recently, the author (reference S; and Imai (reference 4) 
extended the calculations for the circular cylinder "by in- 
cluding the terms involving the fourth power of the Mach 
number. These higher-power terms, neglected in the ear-- 
lier calculations , were f ound to. be of considerable impor- 
tance as the local velocity of sound is apt) reached on the 
surface of the cylinder. It has therefore been thought 
worth while to ext end the calculations ' in a similar manner 
for. the flow past a sphere. 
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^r.§liMn,§TX„il® jelopme The flow is assumed to "be 

uniform at a -^reat distance from the sphere and the mo- 
tion to "be everywhere irrotational and steady. Then, with 



c^-=^^ • CD- 

dp 

the eq^uations of motion reduce to 

. 0^ f = - i ' \ 

.where r is the fluid Telocity; c, the local velocity 
of sound; the pressure; and p, the density. Then, 

assuming the adiahatic relationship "between p and ..p, 
it follows "by integration of equation (2) that 

2 V - 1 p 2 Y - 1 Po : : 

and from . equation (l) that 



where U is the velocity of the undisturhed stream; Pq » 
pQ > ai^d c^ , the corresponding quantities in the undis-*; 
■jtur'bed stream; M (- V/cq), the Mach number; and the 
i^atiio of the specific heats . of the-fluid. ' 

Since the fluid motion is irrotational , there exists 
a velocity potential ^ a.nd the equation "of continuity 
'^B^ "be written as ^ 



-.2 , ^2 , 2 , 

+ + ^ ^ ^ J- 

dJ^ Sz^ 2c^ 



^ll + 111 + ' ^ll ^l) ( 5 ) 
d X 3x d y 3y 3 z dz/ - 



Let r, e, and denote space polar coordinates and sup- 
pose the origin of r to "be at the center of a sphere of 
radius a and the initial line of G to he parallel to 
the direction of the stream. Then, desiijnating "by r the 
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ratio r/a, Tjy f the ratio //Ua and » "by v the ratio 
v/lT and taking into account the fact that the flofrs in 
ail merxdian planes tp = constant are similar, equation 
\^) Decomes ^ 



X - 1 wS 

+ Z M'^* 

where 



1 + ™^ (X - v^) 



Br \^ ary 7^ aH 



(1 - fi^) 1^ 



(7) 



and 



IJ- = cos 0 



. It is now assumed that f can he developed as a power 
series in M"^ (reference 4) so that 



Since 



then 



where 



= To^- + v.^M^ + v/M* + .... 



^ar ar 'r^ 36 36 ^ 



(8) 



(9) 



(9a) 



(9*) 



2 ^ 



Var ar r^ 36 36 / 



1 y3)lS: 



r^ V3T 



» • ■ • 



• • • • t • • • • • ♦ 



(9c) 
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men tjhese expre ssions :f or ^ :-0amM::-:Y^ are inserted into 
equation (6) and the coefflcients^^^^^^^^^^ powers of M 

on tota sides are equated.,. 



• • ♦ • ■ « ■ 



Prom these equations, any <5iven approximation clearly 
depends only on the preceding appi'oxiniations of which the 
ixrst one IS the solution of Laplace equation Ad ^ 0 
for an mcompres sihle fluid. 

-^^^^-^i^-^VapproMmBM^n^ Equation ( 10a) is the dif- 
ferential equation for the r©locl*t^:pot ential for the 

flor of an incompressihle, nonVisooiis fluid and may Ijg 
wri 1 1 en as j 

er \^ Br-J ^ 5]i - ^ > J " ^ 

2?he solution of this equation is known to he 

° n=0 V ^ : j.n+i7 n^i^^ 

where Pj^(p.) is Le?endre»s polynomial cf order n and 
An and are arbitrary constants. In the case of a 

sphere of radius" a, supposed fixe^ in a stream of uni^ 

f±7^ III'' I boundary conditions to he sati^- 

lied are as follows! ■ ' 

¥0 

- dr"" ~ ^o^^al velocity - 0, at the surface of the 

sphere 
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and 

- -— = - cos Q, at infinity 

These conditions limit the form of the solution to 

j^o = (^^"^ + cos 9 

where cos 8 = P^dji), Inserting the "boundary conditions, 

A3^= 1 and B,, = l/2 

and, therefore, 

io = + 2W Px (i^) (11)^ 

SM.-.s£Cond_s:DT2xe3ciffi^j^n . - From equations (ll) and 
(9a), it follows that 

= 1 + + \ r-® + (-3r-^ + | r"^) jx^ 

or since Po(m.) = 1 and (m,) - | ji.^ - "I 



^0 



Then , a sijapXe calculation it i s :found f rom equation 
(10b) that 



18 ^-"^ 



(12) 



where 



\ (|i) = I ^1 



""^"^ " .2 ^ 2 
Now, a particular ..inteijral of the. equation 



r dr A dr/ r"^ dii 



(1 - M^^) 1^1 = r^Pii(M') 
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d = . . ^ILL 

(m + 2) (in + 3) - n(n + ij 



Hence, the soluti<^n is ^iven "by 



2 

n = 0 



ra+j 



'L " (m+2) (m+3) -n(n+l) 

except when ' . 



J 



(13) 



= n - 2 or m = - (n + 3 ) 
Accordingly, the solution of eijuation (r?) is 

" {-To '"'-lo ^""'-^ile -"0^(^> 



S-Lnce already satisfies the necessary boundary condi- 
tions of the pro-blem, the higher approximations, c^. , 
fz* ... must satisfy the conditions 

■ 9r . .' ar.. ' - . .-vi" • 

lation!^ " ^ ''''^ " Hence, after a simple calcu- 

■^1 = As = 0 and 3i ^ -g, = 

I'-hersf ore , ; ^ ■ ' " ■ ^ ' oS^^^^ 

: ^ g r-a _1 +^ r-e) („) 

(- To ^ i - + -I- (,) (14) 

f.?,') -7f-T--ff-f^^^fi-^isa±^^^^^^ Sulstitutin? from equations 
en,!,;/ ^''^ rUht-hana eide of 

iTttlnl^ll''^' " '""""^ straightforward calou- 




V 25 ^' 3-85 ^ 154P^ 1925^ '"15^^' '^ eieT V 

4- f p. r-^ - ?ii ^ §36 ^-o ^ 112 -lo . 3411 12631 ^-13^ 929 

: \ ^ 11 3 77 23 ^ 77 231 ^ ^ 3696 J^o^^^ 



o 



(15) I 

The complste solution of this equation is obtained by means of equation (13) together i?ith the 



Doundary conditions §^=0 for both r = 1 and r = » and is as follows: 
/'a -"^[V 70'' ^24^ 560 2800 r^^^ 

V 10 176 1960 , 6600 ;*3^^^ 

\7 52 1120 8512 

^r-5- J:£r-7+ 121?.r-^+ §5^r-io^ 3049 7789 -lA^. . . 

V 150 1926 3080 84700 * " 18480 3 69600 ' /^^ ^^^ 



~5 

r 




r~3 - ^ 




52 


1120 


49 
150 




243 ^7 . 
1926 


53 
150 




150- 



O 

o 

a 

o 



SI 

rvj 



929 ~ 

■r ■ 



- ^Vf^^., 



55 7150 64680 66100 ■ .&„;/r 

V42 11 3 77 728 1540 , 18480 561792 ^^'50^; 

+ Bir-^Jifvi-) -^Bsr-^PsC^i) +B5r-^Pg(M,) (16)' 
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where 



= 0,03566(7^1)4.0.32614 
B3 = 0.02268(Y-l) + 0.74107 

B5 -0.18261(7.^1) + 0.21216 

calcullle^ll^l't"" easily 
o? ^""^ velocity of the fluid at the surface 

of the sphere is ?iven oy . 

_ 1 3 . n 1 

~ 5^0 = 2 6 ■' 7040 ^^^^ ^^""^ - 1215^ sin 36 

+ (0.10572 sin 8^ 0. 16008 sin 38 + 0. 06434 sin 50 OM^ 

^ CY-DCO.Oliea sin 6 - 0,02475 sin 364-0.02582 gin 56 )m' 

" r ie ^ + Tolo (889 sin 6 - 1215 sin 36 )M^ 

+ (0.11048 sine - 0.17018 sin 36 + 0.07489 sin 5e)M* + .. 

„r,^,^.-~-^T~-----~'^^^^^-^£-^-'S&-M^oh number for thp 

offouAd It''' ^'^^ yelocity of ti^rfluITr^^iTpni velocity 
of sound at any point in the field, a change la the trpe 

wUl r>,»^!?^' potential flow :,=ay no longer exist? It 

U/c • at whf^r;.^' °f interest to determine the value of 
local It . 't" ^^^"'^"y °f the fluid just equals -the 

- .eiooity Cith ^Jui^l-.^' M^r.^n:) ;L%!i%::l; 

'■max = 1.5 + 0.31307M^ + O.^eSSS M* + ... (is) 

polnt^°h;re%'hf'''*^"^-?-^ =* the velocity of sound at a 
oarteiocitv^f velocity of the fluid is equal to the lo- ' 
velocity of sound: that is. at a -ooint where 
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it follows from equation (4) tliat 

2c 2 



+ 1 \ 



Y 1 
2 



In conformity with the usage in this paper, c* 
replaced T?y 5* (« r*) and the equation (19) Tjecomes 



(19) 



It is to "be noted that this equation is essentially 
Bernoulli 's equation and does not depend on the shape of 
the hody Immersod in the fluid. 



is 



2 

Y + 1 



Y + 1 



(20) 



The so-called critical value of U/cq (= kcr 
floff past a sphere is then ohtaihed "by putting t 



) for 



max 



the 



Tallies I and II show, respectively, the values of Vj 
and V* calculated from equations (18) and (20) for vari- 
ous^yalues of the Mach numher; These values are also' shown 
in iigure 1, and the points at which the curves intersect 
give the corresponding values of the. cr,itical Mach numher 
Mcr, The critical values are» respe>cti vely , for the sec- 
ond and the third approximatiohsv U^j. 0.587 and 
0.573. The corresponding critical values of M for 
case of an infinitely long circular cylinder are M^r 
0.420 and Mcr = 0.409 (reference 3). 



^cr 
the 



The_ ve 1 0 c i t y_aM_ t he_p r e^ sur e_M^^ .-If Po , 

Po. and Co are the pressure, the density, and the veloc- 
ity of sound, in the undlsturhed stream, then the density 



P of 
given 



the 
"by 



fluid at a point where the velocity is 



Po 



1 + 



(1 r v^) 



I ■ 
1Y-.I 



(21) 



and the pressure p at: this jboiht is 



or 



Po \PoJ ~ 



^ 1 



(1 - vs) 



Y 
Y-l 



(22) 
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lor an incompre g si ble fluid, .-M = 0 and. expression (23 ) 
reduces to 



P - 



T -g = 1 - v2 (24a) 

2 ° 

and, for a compressi-ble fluid'with Y = 1,408 and M = 
.0*57, . ■ . 



p - Pq . ■ r . s. 

= 4.372 J 1 + 0,06.628(1-^2) . - 'll (24b) 

I Po II J / 

The values of v to l^e used in equations (24) are o"b- 
tained from equation (17 ).. Thus . for the incompre s siljle 
f lui d 

T = I sin 0 (25a) 

and, for the second and the third approximations to the. ; 
compressihle .fluid with M = 0.57.,. 

. V ~ 1.54564 sin 0 - 0.05607. sin 30 . (25h) 
V .= 1.55731 sin 0 « 0.07404 sin 30 . 

■ + 0.d0791 ■ sin 5-0^^ (26c) 

Tahle III lists the values of v computed from these 

formulas and tahle lY -^ives the corresDonding values of 

p — . p , " 
_ — obtained from equations (24)^ Ji£;ures 2 and 3 

2 Po;.U ; ... . . 

show, respectively, the graphs of the velocity and Uhe 
pressure distributions of tables III and IV. 

Lan^ley Memorial Aeronaut! ca-l Laboratory, 

National Advisory Committee for Aeronautics, 
Lan^ley Pield, Ta. , May 4 , 1940. . 
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i 



M 








max 




Second 


■ Tlaird 




approximation 


approximation 


0 


1.5 




' .1 


1.50031 ' 


;: : 1.50317 


.2 


1,51252.. 


1.51309 


.3 


1'.52818 • 


1,53106 • 


.4 


1,55009 


1.65919 


• 5 


1 •57827 


1.60049 


.6 


1,61271 


1,65878 


.7 


1.65S4 


3. . 73 877 



TABLE II 



M 


0,1 


0.2 


0,3 


0.4. 


0.5 


0,6 






.9vl228 


4.5753 


3,0656; 


2.3152 


1,8686 


1.5737 : 


.'■a,>3^'W5;:-' 



13 



•TABLE vliEAl 



Incompressi'ble 
0 

,13074 

,26048 ■ 

,38823 

.51303 

.75000 

.96419 

1*14906 

1*29905 

1,40954 

1,47722 

1.49429 

1^50000 



Second 


Third 


app r 0 X i ma t i o n 


approximation 


0 , 


0 


,12021 


.11991 




; 25946 


.36039 


.35835 


.48008 


.47630 


: •71675 


.70857 


.94496 


.93420 


' !• 15599 


1»14851 


1.33867 


1.34183 


1.48046 


1.49903 


r.-67073 


1*60285 


1.59392 


1.63005 


1.60172 


1.63925 
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Third 

a T) "0 T ft V "l Tlri a'i' A >\ vt 


: 0 


ItOOOOO 


1.08381 


1 • OS 38-1 


5 
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1,06696 
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. 9 321 


1.01620 
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.93301 


.93600 




•73680 

* T*.^ *J 


. 81843 


.82287 


30 


.43750 


.50566 


.51839 


40 


.07035 


. .10801 


.12860 


50 


-.32034 


-.32717 


31098' ■•' 


60 


-.68752 


-.74203 


-.74979 ' 


70 


-.98679 


-1,08069 


-1.12679 


80 


-i«l8216 


-1.30023 


-.1.37912 


85 


-1.23289 


-1 . 35682 


-1.445 79. 


90 


^1*25000 


-1.376 04 


-1.46833 
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Figs. 1,2, 3 
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Figure I.- Critical value of the Mach number for a f^phere 
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Figure..2.- Velocity distribution on Figure 3.- Pressure distribution on 
the surface of a sphere. the surface of a sphere* 
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